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Abstract 

The beta Weibull distribution was introduced by Famoye et al. (2005) and 
studied by these authors. However, they do not give explicit expressions for the 
moments. We now derive explicit closed form expressions for the cumulative distri- 
bution function and for the moments of this distribution. We also give an asymp- 
totic expansion for the moment generating function. Further, we discuss maximum 
likelihood estimation and provide formulae for the elements of the Fisher informa- 
tion matrix. We also demonstrate the usefulness of this distribution on a real data 
set. 

Keywords: Beta Weibull distribution. Fisher information matrix, Maximum likeli- 
hood. Moment, Weibull distribution. 

1 Introduction 

The Weibull distribution is a popular distribution widely used for analyzing lifetime data. 
We work with the beta Weibull (BW) distribution because of the wide applicability of 
the Weibull distribution and the fact that it extends some recent developed distributions. 
This generalization may attract wider application in reliability and biology. We derive 
explicit closed form expressions for the distribution function and for the moments of the 
BW distribution. An application is illustrated to a real data set with the hope that it will 
attract more applications in reliability and biology, as well as in other areas of research. 

The BW distribution stems from the following general class: if G denotes the cu- 
mulative distribution function (cdf) of a random variable then a generalized class of 
distributions can be defined by 

F{x) = lGi.){a,b) (1) 
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for a > and 6 > 0, where 



' ^ B{a,b) B{a,b) 

is the incomplete beta function ratio, By{a, b) is the incomplete beta function, B{a, b) = 
r(a)r(5)/r(a + b) is the beta function and r(.) is the gamma function. This class of 
generalized distributions has been receiving increased attention over the last years, in 
particular after the recent works of Eugene et al. (2002) and Jones (2004). Eugene et al. 
(2002) introduced what is known as the beta normal distribution by taking G{x) in ([T]) 
to be the cdf of the normal distribution with parameters fi and a. The only properties 
of the beta normal distribution known are some first moments derived by Eugene et al. 
(2002) and some more general moment expressions derived by Gupta and Nadarajah 
(2004). More recently, Nadarajah and Kotz (2004) were able to provide closed form 
expressions for the moments, the asymptotic distribution of the extreme order statistics 
and the estimation procedure for the beta Gumbel distribution. Another distribution 
that happens to belong to ([T]) is the log F (or beta logistic) distribution, which has been 
around for over 20 years (Brown et al., 2002), even if it did not originate directly from 

While the transformation ([T]) is not analytically tractable in the general case, the 
formulae related with the BW turn out manageable (as it is shown in the rest of this 
paper), and with the use of modern computer resources with analytic and numerical 
capabilities, may turn into adequate tools comprising the arsenal of applied statisticians. 
The current work represents an advance in the direction traced by Nadarajah and Kotz 
(2006), contrary to their belief that some mathematical properties of the BW distribution 
are not tractable. 

Thus, following ([T]) and replacing G{x) by the cdf of a Weibull distribution with 
parameters c and A, we obtain the cdf of the BW distribution 

F{x) = /i-cxp{-(Ax)<=}(a,&) (2) 

for x>0, a>0, 6>0, c>0 and A > 0. The corresponding probability density function 
(pdf) and the hazard rate function associated with ([2]) are: 

fix) = -g^x^-^exp{-6(Ax)^}[l - exp{-{Xxr}r-\ (3) 

and 

, , cA'^x^^-'^expl— fe(Ax)'^|fl — exp| — (Ax)'^ll''^^ 
r{x) = ' . (4) 

-Dl-cxp{-(A2;)<=}lO, 0) 

respectively. Simulation from is easy: if i? is a random number following a beta 
distribution with parameters a and b then X = {— log(l — B)}^^^/X will follow a BW 
distribution with parameters a, b, c and A. Some mathematical properties of the BW 
distribution are given by Famoye et al. (2005) and Lee et al. (2006). 

Graphical representation of equations and (jl]) for some choices of parameters a 
and b, for fixed c = 3 and A = 1, are given in Figures [T] and [21 respectively. It should 
be noted that a single Weibull distribution for the particular choice of the parameters c 
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and A is here generalized by a family of curves with a variety of shapes, shown in these 
figures. 

The rest of the paper is organized as follows. In Section 2, we obtain some expansions 
for the cdf of the BW distribution, and point out some special cases that have been 
considered in the literature. In Section 3, we derive explicit closed form expressions for the 
moments and present skewness and kurtosis for different parameter values. Section 4 gives 
an expansion for its moment generating function. In Section 5, we discuss the maximum 
likelihood estimation and provide the elements of the Fisher information matrix. In 
Section 6, an application to real data is presented, and finally, in Section 7, we provide 
some conclusions. In the appendix, two identities needed in Section 3 are derived. 




2 Expansions for the distribution function 

The BW distribution is an extended model to analyze more complex data and generalizes 
some recent developed distributions. In particular, the BW distribution contains the 
exponentiated Weibull distribution (for instance, see Mudholkar et al., 1995, Mudholkar 
and Hutson, 1996, Nassar and Eissa, 2003, Nadarajah and Gupta, 2005 and Choudhury, 
2005) as special cases when 6=1. The Weibull distribution (with parameters c and A) is 
clearly a special case for a = b = 1. When a = 1, ([3]) follows a Weibull distribution with 
parameters Xb^^^ and c. The beta exponential distribution (Nadarajah and Kotz, 2006) 
is also a special case for c = 1. 

In what follows, we provide two simple formulae for the cdf ([2]), depending on whether 
the parameter a > is real non-integer or integer, which may be used for further analytical 
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Figure 2: The hazard function (jl]) of the BW distribution, for several values of parameters 
a and b 

or numerical analysis. Starting from the explicit expression for the cdf ([2]) 

^(^) = 7^ exp{-b{Xyr}[l-exp{-{Xyr}r-'dy, 
B{a,b) Jo 

the change of variables {XyY = u yields 

1 |.(Ax)- 



If a > is real non-integer we have 

oo 

r(a-j)j! 



j=o 



It follows that 



5(a, 


&) 


1 








1 





r(a-j)j! 



^^^^ R(n h\ L ^2^ r(n - 

j=0 



oo 



j=0 



oo 



5(«,5)^r(a-j)j!(& + j) 



;-l)Jr(a) _ ^_(fe+j)(Ax)-|_ 
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Finally, we obtain 

For positive real non-integer a, the expansion may be used for further analytical 
and/or numerical studies. For integer a we only need to change the formula used in ([5]) 
to the binomial expansion to give 

When both a and b = n — a + 1 are integers, the relation of the incomplete beta 
function to the binomial expansion gives 



^(^) = E C") [1 - exp{-(Ax)^}]^exp{-(n - j)(Ax)^}. 

j=a ^"'-^ 

It can be found in the Wolfram Functions that for integer a 



(8) 



T( h^ 1 (i-i/)^ f^ r(6 + jV 



and for integer b, 



b-l 



^T{a + j){l-yy 



F a ^ j! 



Then, if a is integer, we have another equivalent form for 



exp{-6(A,T)^} ^ T{b + j 
For integer values of b, we have 



r(6) ^ j! 



^(.) Jl--PW>l° |:I(^exp{-.(A.n. (10) 

Finally, if a = 1/2 and 6 = 1/2, we have 

= — arctany^exp{(Axyy~-T. (11) 

TT 

The particular cases ([HD and (|TUI) were discussed generally by Jones (2004), and the ex- 
pansions (!6l)-( fTTl) reduce to Nadarajah and Kotz's (2006) results for the beta exponential 
distribution by setting c = 1. Clearly, the expansions for the BW density function are 
obtained from ([6]) and ([7]) by simple differentiation. Hence, the BW density function can 
be expressed in a mixture form of WeibuU density functions. 



^http : //functions .wolf reun. co m/GammcLBetaErf /BetaRegularized/03/01/ 
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3 Moments 



Let X be a BW random variable following the density function ([3]). We now derive 
explicit expressions for the moments of X. We now introduce the following notation (for 
any real d and a and b positive) 



POO 

Sd,b,a = / x'^-^exp(-6x){l - exp{-x)Y-^dx. 
Jo 



(12) 



The change of variables x = — log(z) immediately yields Si^b,a = B{a,b). On the other 
hand, the change of variables x = {XyY gives the following relation 

/•oo x -7 

/ y-'-'expi-biXymi - exp{-iXyY}r-'dy = (13) 

Jo c 



from which it follows that 



S'7/,,,,, = E(a,6)A^-^E(X^-'=) 



or, equivalently, for any real r 

relating Sr/c+i,b,a to the rth generalized moment of the beta Weibull. 

First, we consider the integral ( JT2l) when d is an integer. Let f/ be a random variable 
following the Beta(6, a) distribution with pdf /[/(•) and W = — log(?7). Further, let 
F[/(-) and -Fiy(-) be the cdf's of U and W, respectively. It is easy to see that -FV(x) = 
1 — Fu{e^^). Further, by the properties of the Lebesgue-Stiltjes integral, we have 

/oo poo 
x'^-^dFwix) = / x'^-^e"'fuie"=)dx 
-oo Jo 

oo 



B{a,b) 

Sd,b,a 

B{a,b)' 



Thus, the values of Sd,b,a for integer values of d can be found from the moments of W if 
they are known. However, the moment generating function (mgf) Mw{t) = -E'(e*'^) of W 
can be expressed as 



Mw{t) = E{U-') = -/— l\'-'-\l-xr-'dx 

B{a,b) Jo 

B{b-t,a) 



B{a,b) 



This formula is well defined for t < b. However, we are only interested in the limit t ^ 
and therefore this expression can be used for the current purpose. We can write 



Sd,b,a = B{a, b)E{W''-') = B{a, b)Ml}-'^ (0) = ^B{b - t, a) 



(15) 

t=o 
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From equations (fT^ and (ITSl) for any positive integer k we obtain a general formula 

(16) 



E(X^') = J-, ^B(b - t, a) 



t=o 

As particular cases, we can see directly from ( ITSl) that 

S,,b,a = B{a, b), S2,b,a = B{a, b) (a + 6) - ^ (6)} 

and 

53,fe,„ = B (a, 6) [V-' (6) - ^' (a + 6) + {^A (a + 6) - ^ (6)}'], 
and by using (|T^ we find 

^™ ^ {V;(a + &)-^(fe)} 

and 

, ^'(b)-i,'{a + b) + {^{a + b)-^{b)Y 
E{X ) = . 

The same results can also be obtained directly from (fTGll . 

Note that the formula for Si^b,a matches the one just given after equation (12). Since 
the BW distribution for c = 1 reduces to the beta exponential distribution, the above 
formulae for E{X'^) and E{X'^'^) reduce to the corresponding ones obtained by Nadarajah 
and Kotz (2006). 

Our main goal here is to give the rth moment of X for every positive integer r. In 
fact, in what follows we obtain the rth generalized moment for every real r which may 
be used for further theoretical or numerical analysis. To this end, we need to obtain a 
formula for Sci,b,a that holds for every positive real d. In the appendix we show that for 
any d > the following identity holds for positive real non-integer a 

and that when a is positive integer 



j=0 

is satisfied. 

It now follows from (fT7|) and (HM that the rth generalized moment of X for positive 
real non-integer a can be written as 

r(a)r(r/c+l) f. (-1)^- 
^ ^ A^i?(a,6) j-^J(a-j)j\ib + jY/^+^- ^ ' 

When a > is integer, we obtain 

_ r(r/c+l) '^fa-l\ i-iy 
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When a = b = 1, X follows a Weibull distribution and (j20l) becomes 



E{X ) - —^^ , 

which is precisely the rth moment of a Weibull distribution with parameters A and c. 
Equations (fT6l) . (fT9l) and ( |20l) represent the main results of this section, which may serve 
as a starting point for applications for particular cases as well as further research. 




Figure 3: Skewness of the BW distribution as a function of parameter a, for several values 
of parameter b 

Graphical representation of skewness and kurtosis for some choices of parameter b as 
function of parameter a, and for some choices of parameter a as function of parameter 
b, for fixed A = 1 and c = 3, are given in Figures [3] and HI and [5] and El respectively. 
It can be observed from Figures 3 and 4 that the skewness and kurtosis curves cross at 
a = 1, and from Figures 5 and 6 that both skewness and kurtosis are independent of b 
for a = 1. In addition, it should be noted that the Weibull distribution (equivalent to 
BW for a = 6 = 1) is represented by as single point on Figures 3-6. 
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Figure 4: Kurtosis of the BW distribution as a function of parameter a, for several values 
of parameter b 




Figure 5: Skewness of the BW distribution as a function of parameter b, for several values 
of parameter a 
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Figure 6: Kurtosis of the BW distribution as a function of parameter b, for several values 
of parameter a 



4 Moment Generating Function 

We can give an expansion for the mgf of the BW distribution as follows 



M{t) 



B{a,b) Jq 



exp{tx)x^-hxp{-b{Xxy}[l - exp{-{XxY}]''-^dx 



+1 

Y- 



X 



r+c—l 



exp{-6(Ax)'=}[l - ex^{-{Xxy]Y-^dx 



cX' 



B(a, b) ^ r! c 



rlc-\-\,h,ai 



where the last expression comes from (fT3l) . For positive real non- integer a using (fT7|l we 
have 



M{t) 



Via) tT(r/c+l)(-iy 



EE 



{21) 



and for integer a > using ( fTSl) we obtain 



a-l 



a — 1 



(22) 
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Note that the expression for the mgf obtained by Choudhury (2005) is a particular case 
of (pTl) . when a = 9, \ = 1/a and 6=1. When c = 1, we have 



M{t) = y-^-J—-S, 

X 



r=0 ^ ' ^ 



tx—h\x/-\ „—\x\a—l 



1 - e-^^'Y-'dx. 



B{a,b) 

Substituting y = exp(— Ax) in the above integral yields 

and using the definition of the beta function in fl23l) we find 

B{a,b) ' 

which is precisely the expression (3.1) obtained by Nadarajah and Kotz (2006). 



5 Estimation and information matrix 

Let y be a random variable with the BW distribution ([3]). The log-likelihood for a single 
observation y of Y is given by 

i{X,c,a,b) = log(c) + clog(A) + (c -l)log(y)- log {i?(a, 6)} -6(Ay)'= 
+ (a-l)log[l-exp{-(Ay)^}]. 

The corresponding components of the score vector are: 

- {V' (a) - 7/> (a + 6)} + log |l - e-(^^)^| , (24) 



da 



^ = -{^(6)-^(a + 6)}-(Ay)% (25) 



(^ ^ (X ^ , (a-l)(A|/)'^log(Ay)e-(^^)^ 
— = - + log(Ai/)-6(Ai/) log(Ay) + i_e-(\yr ^^^^ 

and 

di c bc,^ c(a-l)(XyYe-^^y^^ 

The maximum likelihood equations derived by equating 0241) - 02 71) to zero can be solved 
numerically for a, b, c and A. We can use iterative techniques such as a Newton- Raphson 
type algorithm to obtain the estimates of these parameters. It may be worth noting from 
E{di/db) = 0, that ([25]) yields 

^(y^c) ^ ^{a + b)-ij{b) 
11 



which agrees with the previous calculations. 

For interval estimation of (a, b, c, A) and hypothesis tests, the Fisher information ma- 
trix is required. For expressing the elements of this matrix it is convenient to introduce 
an extension of the integral ffT^ 

POO 

Td,b,a,e= / x"-ie-^^(l-e-^)""'(logx)^rfx, (28) 
Jo 

so that we have 

Td,b,afi = Sd,b,a- 

As before, let W = — log(t/), where U is a random variable following the Beta(6, a) 
distribution, then 

a;'^-ie^"(l-e-")"-^(logx)"rfa; 



poo 




J —oo 




1 




B{a. 




Td,b, 


a,e 


B{a. 


,by 



Hence, the equation 

Td,b,a,e = B{a, b)E[W''-'{\og{W)}% (29) 

relates Td^b,a,e to expected values. 

To simplify the expressions for some elements of the information matrix, it is useful 
to note the identities 

5'i,b+2,a — '^Si^b+i,a + Si^b,a = B[a + 2, 6) 

and 

bS2,b,a - (a + 26 + l)S2,b+i,a +{a + h + l)52,6+2,a = B{a + 2, 6), 

which can be easily proved. 

Explicit expressions for the elements of the information matrix K, obtained using 
Maple and Mathematica algebraic manipulation software (we have used both for double 
checking the obtained expressions), are given below in terms of the integrals fll2p and 

m-. 



f^a,a = 'ip' {a) -ip' {a + h) 



T2b+ l,a-l,l C l,a-l 



l^afi = -W [a + b) , Ka^c = 57 T-, fia,X 



cB{a,b) ' "'^ XB{a,b) 



1 1 f-L\ in < 7,\ ^2 6 a,l C 5*2 6 a 

Kb,b = V [b)-V {a + b) , Kb^c = —fT-f — ET' = 7- 



cB (a, b) ' ' XB (a, b) 
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+bT2,b,a-2,2 - (a + 2 6 - 1) T2,b+i,a-2,2 + {a + b - l) r2,6+2,a-2,2} , 
+&^2,6,a-2,l - (a + 2 6 - 1) T2,b+i,a_2,l + (o + 6 - 1) T2,6+2,a-2,l} 



A2^A2S (a, 6) 

The integrals ^ and Tij k^i in the information matrix are easily numerically determined 
using MAPLE and MATHEMATICA for any a and b. 

Under conditions that are fulfilled for parameters in the interior of the parameter 
space but not on the boundary, the asymptotic distribution of the maximum likelihood 
estimates d,b,c and A is multivariate normal N^lO, K'^). The estimated multivariate 
normal A^4(0, K^^) distribution can be used to construct approximate confidence intervals 
and confidence regions for the individual parameters and for the hazard rate and survival 
functions. The asymptotic normality is also useful for testing goodness of fit of the 
BW distribution and for comparing this distribution with some of its special sub-models 
using one of the three well-known asymptotically equivalent test statistics - namely, the 
likelihood ratio (LR) statistic, Wald and Rao score statistics. 

We can compute the maximum values of the unrestricted and restricted log-likelihoods 
to construct the LR statistics for testing some sub-models of the BW distribution. For 
example, we may use the LR statistic to check if the fit using the BW distribution is 
statistically "superior" to a fit using the exponentiated WeibuU or WeibuU distributions 
for a given data set. Mudholkar et al. (1995) in their discussion of the classical bus-motor- 
failure data, noted the curious aspect in which the larger EW distribution provides an 
inferior fit as compared to the smaller WeibuU distribution. 

6 Application to real data 

In this section we compare the results of fitting the BW and WeibuU distribution to the 
data set studied by Meeker and Escobar (1998, p. 383), which gives the times of failure 
and running times for a sample of devices from a field-tracking study of a larger system. 
At a certain point in time, 30 units were installed in normal service conditions. Two causes 
of failure were observed for each unit that failed: the failure caused by an accumulation 
of randomly occurring damage from power-line voltage spikes during electric storms and 
failure caused by normal product wear. The times are: 275, 13, 147, 23, 181, 30, 65, 10, 
300, 173, 106, 300, 300, 212, 300, 300, 300, 2, 261, 293, 88, 247, 28, 143, 300, 23, 300, 80, 
245, 266. 



and 
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The maximum likelihood estimates and the maximized log-likelihood Ibw for the BW 
distribution are: 

a = 0.0785, b = 0.0659, c = 7.9355, A = 0.004987 and Ibw = -169.919, 

while the maximum likelihood estimates and the maximized log-likelihood Iw for the 
WeibuU distribution are: 

c = 1.2650, A = 0.005318 and !w = -184.3138. 

The likelihood ratio statistic for testing the hypothesis a = b = 1 (namely, WeibuU 
versus BW distribution) is then w = 28.7896, which indicates that the Weibull distribu- 
tion should be rejected. As an alternative test we use the Wald statistic. The asymptotic 
covariance matrix of the maximum likelihood estimates for the BW distribution, which 
comes from the inverse of the information matrix, is given by 

488130.870 87.9136383 \ 

4009.69885 -135.603333 

58517447.8 -16222.8149 ' 
16222.8149 6.19530131 / 

The resulting Wald statistic is found to he W = 38.4498, again signahzing that the BW 
distribution conform to the above data. In Figure 7 we display the pdf of both Weibull 
and BW distributions fitted and the data set, where it is seen that the BW model captures 
the aparent bimodality of the data. 



/ 8699.35364 
4743.69977 
-488130.870 
\ 87.9136383 



4743.69977 
13079.4394 
-4009.69885 
-135.603333 




Figure 7: The probability density function ([3]) of the fitted BW and Weibull distributions 
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7 Conclusion 

The Weibull distribution, having exponential and Rayleigh as special cases, is a very pop- 
ular distribution for modeling lifetime data and for modeling phenomenon with mono- 
tone failure rates. In fact, the BW distribution represents a generalization of several 
distributions previously considered in the literature such as the exponentiated Weibull 
distribution (Mudholkar et al., 1995, Mudholkar and Hutson, 1996, Nassar and Eissa, 
2003, Nadarajah and Gupta, 2005 and Choudhury, 2005) obtained when 6=1. The 
Weibull distribution (with parameters c and A) is also another particular case for a = 1 
and 6=1. When a = 1, the BW distribution reduces to a Weibull distribution with 
parameters A6^/'^ and c. The beta exponential distribution is also an important special 
case for c = 1. 

The BW distribution provides a rather general and flexible framework for statistical 
analysis. It unifies several previously proposed families of distributions, therefore yielding 
a general overview of these families for theoretical studies, and it also provides a rather 
flexible mechanism for fitting a wide spectrum of real world data sets. 

We derive explicit expressions for the moments of the BW distribution, including an 
expansion for the moment generating function. These expressions are manageable and 
with the use of modern computer resources with analytic and numerical capabilities, may 
turn into adequate tools comprising the arsenal of applied statisticians. We discuss the 
estimation procedure by maximum likelihood and derive the information matrix. Finally, 
we demonstrate an application to real data. 

Appendix 

In what follows, we derive the identities ( |T7I) and ( |T8i) . We start from 



exp(-6x)(l - e ^■) 



which yields 




and substituting z = e ^ gives 




(30) 



For real non-integer a, we have 




oo 

f{x)dx = r(a) ^ 

j=0 



r(a 



(-1) 



Also, for real p > —1 and real q, we have 




(31) 
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Hence, 

and, finally, we arrive at 



/•oo 

/ x'^-^exp{-bx){l - e-'^f-^dx = T{a)T (d) V 
Jo 



which represents the identity (1171) . 

Now, let a > be an integer; then, from (!30|) . we have 



r x'''^f{x)dx = 5^ f"" . ^V-iy |logz|''-^/+^-Vz. 



Using (pTl) we obtain 

a-1 



a - 1\ (-l)^r(rf) 



J / + 



d 



j=0 

and therefore we arrive at 

a-1 



a-l\ (-IV 



which represents the identity ( |T8l) . 

References 

[1] Brown, B. W., Spears, F. M. and Levy, L. B. (2002). The log F: a distribution for 
all seasons. Comput. Statist., 17, 47-58. 

[2] Choudhury, A., 2005. A simple derivation of moments of the exponentiated WeibuU 
distribution. Metrika, 62, 17-22. 

[3] Eugene, N., Lee, C. and Famoye, F., 2002. Beta- normal distribution and its appli- 
cations. Commun. Statist. - Theory and Methods, 31, 497-512. 

[4] Famoye, F., Lee, C. and Olumolade, O., 2005. The Beta- WeibuU Distribution. J. 

Statistical Theory and Applications, 4, 121-136. 

[5] Gupta, R. D. and Kundu, D., 2001. Exponentiated Exponential Family: An Alter- 
native to Gamma and Weibull Distributions. Biometrical Journal, 43, 117-130. 

[6] Gupta, A. K. and Nadarajah, S., 2004. On the moments of the beta normal distri- 
bution. Commun. Statist. - Theory and Methods, 33, 1-13. 



16 



[7] Jones, M. C, 2004. Families of distributions arising from distributions of order statis- 
tics. Test, 13, 1-43. 

[8] Lawless, J. F., 1982. Statistical Models and Methods for Lifetime Data. John Wiley, 
New York. 

[9] Lcc, C, Famoyc. F. and Olumoladc, O., 2007. Beta-WeibuU Distribution: Some 
Properties and Applications to Censored Data. J. Modern Applied Statistical Meth- 
ods, 6, 173-186. 

[10] Linhart, H. and Zucchini, W., 1986. Model Selection. John Wiley, New York. 

[11] Meeker, W. Q. and Escobar, L. A., 1998. Statistical Methods for Reliability Data. 
John Wiley, New York. 

[12] Mudholkar, G. S., Srivastava, D. K. and Freimer, M., 1995. The exponentiated 
WeibuU family. Technometrics, 37, 436-45. 

[13] Mudholkar, G. S. and Hutson, A. D., 1996. The exponentiated WeibuU family: some 
properties and a flood data application. Commun. Statist. - Theory and Methods, 
25, 3059-3083. 

[14] Nadarajah, S. and Gupta, A. K., 2005. On the moments of the exponentiated WeibuU 
distribution. Commun. Statist. - Theory and Methods, 34, 253-256. 

[15] Nadarajah, S. and Kotz, S., 2004. The beta Gumbel distribution. Math. Probab. 
Eng., 10, 323-332. 

[16] Nadarajah, S. and Kotz, S., 2006. The beta exponential distribution. Reliability 
Engineering and System Safety, 91, 689-697. 

[17] Nassar, M. M. and Eissa, F. H., 2003. On the exponentiated WeibuU distribution. 
Commun. Statist. - Theory and Methods, 32, 1317-1336. 



17 



